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Abstract. We present millimeter imaging observations in the 1200 /im continuum of the disk around /3Pic. With 
the 25" beam, the /3Pic disk is unresolved perpendicularly to the disk plane (< 10"), but slightly resolved in 
the northeast-southwest direction (26"). Peak emission is observed at the stellar position. A secondary maximum 
is found 1000 AU along the disk plane in the southwest, which does not positionally coincide with a similar 
feature reported earlier at 850 /im. Arguments are presented which could be seen in support of the reality of 
these features. The observed submm/mm emission is consistent with thermal emission from dust grains, which 
are significantly larger than those generally found in the interstellar medium, including mm-size particles, and 
thus more reminiscent of the dust observed in protostellar disks. Modelling the observed scattered light in the 
visible and the emission in the submm/mm provides evidence for the particles dominating the scattering in the 
visible/NIR and those primarily responsible for the thermal emission at longer wavelengths belonging to different 
populations. 
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1. Introduction 



Since its discovery by IRAS (e.g., IAumann~ 984). the 
/3Pictoris system has presented the prime example of 
a dusty disk around a main sequence star, partly be- 
cause of its high degree of 'dustiness' (Ltr/£* = 
2.5 x 10~ 3 , e.g. |Lagrange et al. 2000| ) and partly because 
of its relatively close distance to the Earth (19.3 pc, 
Crif o et al. 1997(1 . which makes it possible to obtain high 
quality data over the entire spectral range. Recent pa- 
pers reviewing the physics of the disk around (3 Pic 
incl ude those of |Artymowicz 2000| |Lagrange et al. 2000 
and IZuckerman 200ll Considerable uncertainty existed 
regarding the age of the system, but most recent esti- 
mates place the stellar age close to only ten million years 
(12+^Mvr. [Zuckerman et al. 2001jl . This could open up 
the possibility that planet formation (nearing its final 
phases?) might actually become observable. 



circumstellar matter - planetary systems: formation - protoplanetary 



Since these reviews were written, new relevant infor- 
mation has been added to our knowledge of the /3Pic 
system: IQlofsson et al. 200T1 reported the discovery of 
widespread atomic gas in the disk, a result recently con- 
firmed and extended by Brandeker et al. 20021 These ob- 
servations revealed the sense of disk rotation and that the 
northeast (NE) part of the gaseous disk is extending to 
the limit of the observations by Brandeker et ah, viz. to 
at least 17" (330 AU) from the star. Several difficulties 
were encountered with these discoveries, such as the ob- 
served fact that the gas stays on (quasi-)Keplerian orbits, 
although radiation pressure forces in the resonance lines 
should accelerate the gas to high velocities and remove it 
on time scales comparable to the orbital period. 



Send offprint requests to: R. Liseau 

* Based on observations collected with the Swedish ESO 
Submillimeter Telescope, SEST, in La Silla, Chile. 



This needs to be addressed in the context of the ori- 
gin and evolution of the gas and dust, i.e. whether one or 
both components are presently produced in situ in the 
disk or whether they (at least to some degree) consti- 
tute 'left-overs' from the star formation process, being of 
primordial origin. Tak euchi fc Artymowicz 2001| have con- 
sidered the interaction of gas and dust in a circumstellar 
disk, the dynamical evolution of which is critically depen- 
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Fig. 1. The normalised SIMBA 1.2 mm image with 8" pix- 
els of the flux calibrator Uranus, the size of which was 3"5 
in diameter at the time of our observations and, hence, ap- 
peared point-like to SIMBA. Offsets in Right Ascension 
and Declination are in arcsec and the derived circular 
Gaussian beam of 25" FWHM is shown in the lower left 
corner. During the scanning alt-az observations, the im- 
age rotates which would smear out any low-level features. 
Contour levels as in Fig.|3 
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Fig. 2. The f3 Pic disk imaged at 1.2 mm with a pixel sam- 
pling of 8". The lowest (dashed) contour corresponds to 
2a and increments are by ler = 3mJy/beam. The center 
coordinates (0, 0) refer to the stellar position and offsets 
are in arcsec. Features discussed in the text are marked 
and the optical disk midplane is shown by the straight 
line at position angle 31°5. The circular Gaussian beam of 
25" FWHM (~ 500 AU) is shown to scale in the lower left 
corner. 



dent on the relative abundance of these species (see also 
|Lecavelier des Etangs et al. 1998| |. Possible observational 
consequences, even relatively far from the central star, 
may become assessible with modern mm/submm cameras. 
With the aim to compare and to extend the results ob- 
tained at 850 /im with SCUBA by Holland et al. 1998, 
we performed imaging observations at longer wavelengths 
and, in this paper, we present the image of the /?Pic disk 
at 1200 /im (1.2mm). 

The SIMBA observations and the reduction of the 
data are presented in Sect. 2. Our basic result, i.e. the 
1.2 mm image of /3Pic and its circumstellar disk, is found 
in Sect. 3, and in Sect. 4 the possible implications of these 
observations are discussed, where also other data are con- 
sulted. Finally, in Sect. 5, we briefly summarise our main 
conclusions from this work. 

2. SIMBA observations and data reductions 

The observations were performed with the SEST Imaging 
Bolometer Array (SIMBA) at the 15 m Swedish-ESO 
Submillimeter Telescope (SEST), La Silla, Chile, dur- 
ing the periods November 17 - 23 and 27 - 29, 2001. 
SIMBA has been developed by the Max-Planck-Institut 
fur Radioastronomie, Bonn, in collaboration with the 
Astronomisches Institut der Ruhr-Universitat Bochum. 
The 37 liquid helium cooled semiconductor elements are 
n-doped silicon chips, mounted on a saphire substrate. 



Projected on the sky, the 37 horn antennae have each an 
HPBW of 23" and are spaced by 44" in a hexagonal ar- 
rangement, covering about 4'. The spectral band pass is 
centered at 250 GHz (1200 /im) and has a full width at 
half maximum (FWHM) of 90 GHz. 

As map centre, the equatorial coordinates of (3 Pic were 
used, viz. 05 h 47 m 17 s l, -51°03'59" (J 2000). Generally, the 
mapping was done in fast mode by scanning 600" in az- 
imuth and 400" in elevation at the rate of 80" per second 
and with a step size in elevation of 8", oversampling the 
beam by about a factor of three. The zenith opacities were 
obtained by means of frequent sky dips and were on the 
average 0.22 at the beginning of our observing run, but 
improved to about 0.15 after a couple of days. The point- 
ing of the telescope was regularly checked using the extra- 
galactic radio source 0537 — 441 and/or the planet Saturn 
and was found to be within a third of the beamwidth. 

Flux calibrations were based on observations of the 
planet Uranus. The brightness temperature of 94.8 K at 
250 GHz implies a Uranus flux density of 38Jy/beam, 
with an uncertainty of 5% (Grim n~&: Orton 1993|l . As is 
evident from Fig.Q] the spatial flux distribution of this 
point source (diameter 3"5) is consistent with a circular 
Gaussian telescope beam, with the measured FWHM( | 
, J_) = (24"9 ± 0"3, 24"6 ± 0"3). The deviation from a cir- 
cularly symmetric Gaussian beam pattern occurs at the 
1% level of the peak value (— 20 db). 
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In total, we recorded 171 maps with 21 hours of in- 
tegration. These data were reduced by making use of the 
MOPSI reduction software package 1 developed by Robert 
Zylka. This involved correcting for atmospheric extinction, 
cosmic rays and the variable sky background, as well as 
producing maps from the fast scanning mode. The sky 
noise was greatly reduced by correlating and removing the 
simultaneous flux level variations of the different bolome- 
ter channels. 

3. Results 

The 1200 /xm image of the /3Pic region is displayed in 
Fig. [21 where the displayed contours are chosen in compli- 
ance with Holland et al. 19981 The flux maximum, desig- 
nated as feature A in the figure, is centered on the position 
of the star (3 Pic (0", 0")- At position angle 237° relative to 
A, an elongation of the emission in the NE-SW direction 
is discernable. For an assumed Gaussian source flux distri- 
bution, the deconvolved minor and major axes are < 10" 
and 26", respectively, the latter corresponding to 500 AU. 
The disk emission extends to at least 55" (1050 AU) in the 
SW direction. In addition, faint emission protrudes south 
toward a 'blob' at (-27", -44"). This feature C, at posi- 
tion angle 211°5, is 52" (1000 AU) distant from the star 
and is thus not (whether real or not) positionally coinci- 
dent with the 850/xm-SW blob of Holland et al. 19981 at 
(-21", -26") and identified as B in Fig.H but we note the 
following important facts: (1) this faint emission, extend- 
ing straight south, is also apparent in the 850 /im image, 
(2) as is evident in Fig.El the position angle to blob C co- 
incides with that of the midplane of the (3 Pic disk inferred 
from optical data (= 31°5 + 180°; IKalas & Jewitt 19951 
|Heap et al. 2000| ) , and (3) faint dust scattered light in this 
direction has been observed far from the star (1450 AU, 
ILarwood fc K alas 2001). Toward the NE, the scattering 
disk has been claimed to extend even further from the 
star, to 1835 AU. 

An about 2a feature at (+48", +32") and pa = 57° 
(1100 AU) is discernable in our image, but which would 
again be significantly offset from a similar blob at (+28", 
+25") in the 850 /xm data. In order to assess the reality 
of these faint features we have divided the raw data into 
different portions and then applied the same reduction 
procedures to these subsamplcs. The result of this exer- 
cise indicates that the (+48", +32") feature is an artefact 
introduced by the noise, as it is not seen in all frames, 
whereas blob C is persistently present in the final sub- 
maps. 

The flux densities of both the (3 Pic disk (^4) and the 
SW features (B and C) are presented in Tabled The result 
for the peak flux obtained by IChini et al. 1991? with a 
single element bolometer at the SEST at 1300 /im, viz. 
24.9 + 2.6 mjy/beam, is in agreement with our array value 

1 http://www.ls.eso.org/lasilla/Telescopes/SEST/SEST.html 

2 We assume that these data refer to (3 Pic and not to a Pic 
as written in their Table 1. 



Table 1. SIMBA 1200 /xm flux densities of the /3Pic disk 



Feature 


Relative 
Offset (") 


F„(1200 fim) 
(mjy/beam) 


Remarks 


A 

A 


In n\ 
(U, U) 


z4.o it o.U 

35.9 ±9.7 


p Jr ic aiSK 

integrated over a radius of 40" 


B 


(-21, -26) 




SW blob ^Holland et al. 1998fl: 

contaminated by j3 Pic disk 


C 


(-27, -44) 


9.7 ±3.0 


SW blob (this paper) 



of 24.3 ± 3.0mJy/beam at 1200/mi (cf. Fig.|3J). The ratio 
of the peak flux A to that of blob C is 2.5 ± 0.8. 

4. Discussion 

4.1. The /3Pic disk 

For widely adopted parameters, the stellar disk subtends 
an angle in the sky of less than 1 mas and the photosphere 
generates a flux density at the Earth of less than 1 mJy 
at 1200 /im. The stellar contribution to our SEST mea- 
surements can therefore be safely ignored. Also, any line 
emission in this band pass is likely to be totally negli- 
gible ( |Liseau fc Artymowicz 1998| ILiseau 19991 and refer- 
ences therein). 

Assuming an opacity law of the form k v = kq (y/vo)P, 
the flux density at long wavelenghts from an optically thin 
source of a certain dust population can be expressed as 

Fu = ^f^+0 /r dust (z)dz . (i) 

C Z Vq J 

If the particles dominating the emission at 850 /im 
and 1200 /im, respectively, can be assumed to yield the 
same integral, e.g. because they share the common tem- 
perature Td U st and/or occupy similar locations in space 
along the line of sight z, the average spectral index (3 can 
be obtained from our 1200 /im and the 850 /xm fluxes by 
Holl and et al. 19981 using 

dlogF„ log (^850/-Fl200) , v 

1 dlog^ log (1200/850) ' 1 ' 

yielding j3 — 0.5 for a point source at the stellar lo- 
cation. The index becomes j3 = 1, if we use the fluxes 
for the extended source, viz. integrated over a radius of 
40" centered on the star. This includes blob B, contribut- 
ing some 20% to the 850 /xm flux. A 'correction' for this 
would again indicate a lower (3 value and we conclude 
that the dust in the /3Pic disk has a shallow opacity 
index, perhaps even below unity (Drill et al. 20001 sug- 
gest (3 = 0.8). This is illustrated in Fig.|3 which dis- 
plays the long wavelength spectral energy distribution of 
(3 Pic, together with weighted Rayleigh-Jeans spectra for 
(3 — 0, land 2, and could mean that the grains dom- 
inating the millimeter-wave emission are different from 
those scattering most efficiently in the visual and the 
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near infrared. This value of (3 is significantly lower than 
those found in the interstellar medium (ISM), where typ- 
ically (3 ~ 2 (Hilde Tjrand 1983(1 . but it is similar to 
that found in protostellar disks (e.g.. lBeckwith et al. 1 990 
|Mannings fc Emerson 1994| |Dutrey et al. 1 996 ) , indicat- 
ing significant differences between the dust particles in 
the ISM and those in the (3 Pic disk. As was also already 
concluded by IChini et al. 19911 the presence of relatively 
large grains is suggested, with maximum radii in excess 
of 1mm (maxa/A<; 1). Intriguing, however, is the exis- 
tence of such large grains possibly as far away as 1000 AU 
from the star (see, e.g., |Takeuchi fc Artymowicz 2001| 
|Lecavelier des Etangs et al. 1998| |. 

4.2. Weak dust features in the 1.2 mm image 

The asymmetric flux distribution displayed in Fig.0 
may be surprising. In agreement with our observations, 
IChini et al. 19911 and IDent et al. 2 000 too were unable 
to detect any emission in the northeast part of the 
disk, where we place a 3er-upper limit on the mass 
of 0.2 M e (see below; the dust temperature at 500 AU 
T dust = 45K, when T dust (r ) = 110 (r/26 AU)" 3 K, see 
|Liseau fc Artymowicz 1998) 1. 

Given the low signal-to-noise ratio (S/N), the reality 
of this lopsidedness is difficult to assess, but asymmetries 
in the /3Pic disk have been noticed also at other wave- 
lengths. For instance, in scattered light, the receding NE 
side of the disk extends much further and is much brighter 
than the SE disk. In contrast, the shorter, approach- 
ing SE disk seems much thicker l|Kalas fc .Tewitt 1 995 
ILarwood fc K alas 2001} . The situation is reversed in 
the thermal infrared (albeit on smaller spatial scales), 
where the SW disk appears significantly brighter and 
more extended than the NE side (|Lagag_ e fc Pantin 1 994 
IWahhaj et al. 20021 |Weinberger et al. 2002| ). This could 
be due to a 'Janus-effect', i.e. the NE being dominated 
by 'bright' dust particles (high albedo, silicates), whereas 
in the SE, the majority of dust grains is 'dark' (high ab- 
sorptivity, carbonaceous?). What would accomplish such 
uneven distribution in the disk is not clear, but large dif- 
ferences in albedo, by more than one order of magnitude, 
are not uncommon, for instance, in solar system material. 
Also, in order to understand the nature of feature C (and 
B) velocity information would be valuable. 

Blob C would be situated in the disk midplane and 
on the second contour in the scattered light image of 
ILarwood fc Kalas 20011 (22 < R < 25 mag/arcsec 2 , see 
also Ka las <^ .Tewitt T9 95 ). No obvious distinct feature 
is seen at its position. However, to be detectable with 
SIMBA at the SEST, any point source at mm- wavelengths 
would not be point-like at visual wavelengths, and its op- 
tical surface brightness could be very low. 

To gain some quantitative insight we ran numer- 
ical models, exploiting Mie-theory, for a variety of 
plausible dust mixtures regarding the chemical com- 
position and grain size distributions (for details, see 




10 100 1000 

Wavelength (/U.m) 

Fig. 3. The SED of the /3Pic di sk from 10 to 1300 /urn. 
Data up to 200 ^m are from fHeinrichs en et al. 1 999 
(52" to 180"). The da tum at 800 /im is from 
Zuc kerman fc Becklin 19931 (integrated over a radius = 
25"). The 850 /Ltm data are from IHolland et al. 19981 and 
refer to a 14" beam and to the integrated flux over a radius 
= 40". The latter should be comparable to our 1200 fim 
point for 40"; the lower is for a 25" beam, and that at 
1300 fim (24" beam) is from IChini et al. 19911 For refer- 
ence, the spectral slopes longward of 100 /im for three val- 
ues of the dust emissivity index /3, discussed in the text, 
are shown by the straight lines. 

IPantin et al. 199"7j) . The equilibrium temperatures were 
found to be in the interval 16 K to 58 K and depend- 
ing on the dust albedo and scattering phase function, 
the predicted integrated scattered light, to be consistent 
with our SEST observations, spans 8 magnitudes. The two 
most extreme cases considered were (1) bright cold dust 
(albedo = 1 at visual wavelengths, 7d ust = 16 K) which is 
scattering isotropically and has albedo = 0.2 at thermal 
wavelengths, which results in a spatially integrated R- 
magnitude of 17.3 and (2) dark warm dust (albedo = 0.02 
in the visual, Tdust = 58 K) which gives rise to 'comet 
scattering', i.e. 14% of isotropic at 90°, and has zero 
albedo at thermal wavelengths, resulting in an integrated 
i?-magnitude of 25.6. These extreme cases are felt to be 
either too optimistic or too conservative and an interme- 
diate case might be more appropriate. Our adopted model 
includes isotropically scattering dust at Tdust = 25 K with 
albedo = 0.2 at visual and zero albedo at thermal wave- 
lengths, yielding an integrated i?-magnitude of 20. An 
about 10" source (R = 25 mag/arcsec 2 ) would thus be 
consistent with both the optical data and our SIMBA 
measurement and a deep i?-band search might become 
successful. Similarly, the integrated 850 /im flux density 
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is predicted to be slightly less than 18mJy and should 
become readily detectable. However, blob C is situated 
outside the figure of Holla ncTet al. 1 998 

Feature C is perfectly aligned with the optical disk 
plane, and its (hypothetical) mass can be estimated from 



dust 



Ky i^(Tdust) 



(3) 



where the distance D is assumed to be that of /3Pic 
and the adopted dust absorption coefficient K250GHZ ~ 
lcm 2 g _1 . This estimate of k v is probably correct within 
a factor of three (Beck with et al. 20001 and references 
therein). For the dust temperature Tdust = 25 K, the dust 
mass is of the order of ten lunar masses (0.16 M^), which 
would be comparable to the dust mass of the disk proper, 
being overall much warmer. Another factor of about three 
uncertainty stems thus from the dust temperature, pro- 
vided 10 K < T dust < 60 K. 

The famous 'SW-blob' in the 850 /im image has re- 
ceived considerable interest by the debris disk community. 
According to IDent et al. 2000| this feature, labelled B in 
Fig-12 is rea l- It is not readily apparent in our 1200 /im 
image, however, presumably due to the combination of 
low contrast and reduced angular resolution. To test this 
idea, we performed numerical experiments, i.e. convolving 
two point sources, at the appropriate positions of A and B 
and with varying flux ratios, with the SIMBA beam. Flux 
ratios, normalised to the peak value, which are consistent 
with our observations are in the range 0.25 — 0.45, with the 
most compelling being about 0.3 (comparable to that for 
blob C). In combination with feature C, blob B accounts 
for the southward bridge seen in Fig. [21 

For radiation mechanisms generating power law spec- 
tra and/or thermal dust emission from blob B, the spec- 
tral slope is given by a = [3 a + 2 — Alogi?A/AlogA, 
where (3a, as before, refers to A — (0", 0"), i.e. the 
/3Pic disk, and where R x = [F V (A) / F V {B)] X and A = 
850 /im or 1300 /im. Because of the relatively low S/N of 
the SCUBA and SIMBA data, the actual flux ratios are 
highly uncertain and, furthermore, calibration uncertain- 
ties and telescope beam effects could potentially introduce 
large errors. The combined observations of blob B suggest 
log (i? 850 /i?i2oo) ~ 0, yielding [3 B ~ (3a, i-e. consistent 
with the spectrum of A. 

5. Conclusions 

Based on 1200 /im imaging observations of the circumstel- 
lar dust of (3 Pic we conclude the following: 

• At 25" resolution, the 1200 /im image of the /3Pic disk 
is slightly resolved in the NE-SW direction, but re- 
mains unresolved perpendicularly (NW-SE). In addi- 
tion, the emission appears asymmetric and extends fur- 
ther SW-S to more than 1000 AU away from the star. 
Maximum emission {— 24mJy) is observed toward the 
position of the star. 



• Combining our 1200 /im map with that by 
Holl and et al. 19981 at 850 /im we infer that the 
dust size distribution in the (3 Pic disk is significantly 
different from that in the general interstellar medium 
and appears more reminiscent of that found in pro- 
tostellar disks. We argue that the thermal emission 
is dominated by dark big particles and that these 
grains constitute a population different from that 
dominating the scattering in the visible/NIR part of 
the spectrum. 

• From the examination and numerical simulations of 
the available data (optical and submm/mm) for the 
southwestern blobs/extensions seen in the SIMBA 
maps we conclude that their reality can at present nei- 
ther be excluded nor can, on the basis of the available 
evidence, their existence be fully confirmed. 
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